arXiv:1506.01470vl [physics.flu-dyn] 4Jun2015 


1 


Comment on “The velocity field due to an 
oscillating plate in an Oldroyd-B fluid” by C.C. 
Hopkins and J.R. de Bruyn [Can. J. Phys. 92, 533 
( 2014 )] 

Ivan C. Christov 


Abstract: We correct certain errors and ambiguities in the recent pedagogical article by Hopkins and de Bruyn. The 
early-time asymptotics of the solution to the transient version of Stokes’ second problem for an Oldroyd-B fluid in a 
half-space is presented, as an Appendix, to complement the late-time asymptotics given by Hopkins and de Bruyn. 

PACS No.; 47.50.d 


Resume: En Fran§ais... 

In a recent pedagogical article, Hopkins and de Bruyn [1] 
present a treatment of the unsteady boundary layer created by 
the oscillatory motion of a plate in a certain viscoelastic fluid. 
Unfortunately, several errors and ambiguities appear in their 
discussion. Given that the goal of [1] is to be an accessible 
and, more importantly, useful introduction to the topic for un¬ 
dergraduate students, who are unfamiliar with the subtleties of 
non-Newtonian fluid mechanics, it is crucial to correct the er¬ 
rors and clarify the ambiguities in [ 1 ]. 

1. In the Introduction and in Section 3 of [1], the authors 
say that they are solving for a “steady-state velocity field.” 
This statement is ambiguous since the final velocity held 
quite clearly depends on time, hence it is unsteady. It is 
more accurate to say that the solution sought and found 
is the late-time (post-transient) velocity held. 

2. The authors mention, in passing, above Eq. (1) in [1] that 
r represents the shear (deviatoric) part of stress tensor 
<j. This is important to emphasize since non-Newtonian 
constitutive relations are, by convention, written in terms 
of T, while it is cr that enters the conservation of linear 
momentum equation. The difference between cr and t 
is the isotropic stress (constitutively indeterminate for 
incompressible fluids), usually identified as the pressure 
p, i.e., cr = —pI-\-T, where I is the identity tensor. Note 
that this is not just a pedantic remark since the general 
form of the right-hand side of Eq. (4) in [1] is V • cr, 
which then makes clear to the reader the origin of term 
-Vp. 

3. Equation (3) in [1] is missing a term. Since, at this point 
in the discussion, the authors have not restricted to in¬ 
compressible fluids or isochoric flows, then the correct 
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form of the upper-convected invariant time derivative of 
a tensor A (see Sec. 3(a) in [2] or Eqs. (8.32.2) and 
(8.32.3) in [3]) is 

l=^-f(M-V)A-(VM)TA-AVM-h(V-M)A, (1) 
at 

where the last term in the equation above is missing from 
Eq. (3) in [1]. 

4. Section 2 in [1] concludes with the following incorrect 
statement; “As an aside, we note that the Cattaneo equa¬ 
tion for heat transport has exactly the same form as ( 2 ) 
[21]. In this case, there are no physical restrictions on 
the values of Xi and A 2 .” In the latter, [21] is [4] below. 
Referring to Eq. (1.2) (“Cattaneo’s equation”) in [4], it 
should be clear that the formal analogy between the con¬ 
stitutive equation relating the deviatoric stress t to the 
rate of strain 7 and the constitutive equation relating the 
heat flux q to the temperature gradient VT only holds 
for the case when A 2 = 0 (using the notation of Eq. (2) 
in [1]). It is also worth noting that the model attributed 
to Cattaneo [5] was suggested earlier by Maxwell [ 6 ]. 
Einally, the value of Ai is not unrestricted in the heat 
conduction context. It is known that if X\ is “too large,” 
then the Maxwell-Cattaneo heat conduction model can 
lead to the violation of the second law of thermodynam¬ 
ics [7]. 

5. Throughout the article, Eq. (4) in [ 1 ] is referred to as “the 
Navier-Stokes equation” but this only true for the New¬ 
tonian constitutive relation. Equation (4) in [1] is a much 
more general relation known as Cauchy’s First Law of 
Motion, expressing the conservation of linear momen¬ 
tum for a continuum (see, e.g., Section III .6 in [ 8 ]). Mean¬ 
while, “the Navier-Stokes equation” usually refers to the 
equation resulting from substituting the Newtonian con¬ 
stitutive relation (the so-called “linearly viscous fluid”) 
into Cauchy’s equation (see, e.g., Eq. (IV.4-25) in [ 8 ] 
and the discussion thereafter). 


unknown 99: 1-3 (2015) 


DOI: m.tmiZxx-xxx 


(c) 2015 NRC Canada 






2 


unknown Vol. 99, 201 5 


6 . Below Eq. (12) in [1], the authors emphasize that the 
“constitutive relation is nonlinear”. While, generally speak¬ 
ing, this statement has some merit, nonlinear terms such 
as u-Vu vanish identically for the assumed unidirec¬ 
tional velocity field given in Eq. (11) of [1]. Moreover, 
the extra assumptions regarding the amplitude of the os¬ 
cillations and the (non-)excitation of higher-harmonics 
are neither used nor needed to obtain the linear ordinary 
differential equation (Eq. (32) in [1]) governing the spa¬ 
tial structure of the post-transient velocity field. Hence, 
we find these several sentences below Eq. (12) in [1] to 
be unclear and potentially confusing for the uninitiated 
reader. 

7. Above Eq. (22) in [1], “the constitutive relation, (3)” 
should be replaced by “the constitutive relation, (15)”. 
Eurthermore, in the top-left entry of the matrix on the 
left-hand side of Eq. (22) in [1], “Tyx — Txy” should be 
replaced by “Tyx -f Txy”. Similarly, “Tyx — Txy” should be 
replaced by “Tyx -f Txy” in the third term on the left-hand 
side of Eq. (23) in [1]. 

8 . In Eigs. 2 and 3 in [1], (dimensionless) values of Xi and 
A 2 up to 10 are used. Although this is fine for illustrating 
the qualitative nature of the curves plotted, it is impor¬ 
tant to note that the two relaxation times are expected to 
be small (especially compared to the inverse plate fre¬ 
quency 1/(0 used as the time scale in [1]), based on 
available experimental data [9]. 

9. Last, but not least, a number of references cited in [1] 
are mathematically erroneous, hence we strongly advise 
readers to avoid these papers. It behooves the commu¬ 
nity to cease the promulgation of these incorrect results 
that plague the modern fluid mechanics literature. Es¬ 
sentially, the mistake can be boiled down to a sloppy en¬ 
forcement of the start-up condition for the transient, un¬ 
steady problem. Eor example, in [10], it was shown that 
Refs. [7,8,14] from [1] contain mathematical errors; in 
[11], it was shown that Ref. [9] from [1] contains math¬ 
ematical errors. Meanwhile, Ref. [15] in [1] makes use, 
without attribution, of the mathematical technique intro¬ 
duced in [12]. Although Refs. [11,13] from [1] do not 
immediately fall within the classes of mathematical er¬ 
rors described in [10, 11, 12, 13], we caution the reader 
against assuming that the results therein are correct. 
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Appendix A: Early-Time Asymptotics 

In this Appendix, we seek a solution to Eqs. (14)-(16)in [1], 
without making the assumption that the velocity is separable, 
namely we do not assume u{y,t)= f{y)e^' as in Eq. (12) in [1]. 
The exact solution to the transient version of Stokes’ second 
problem for an Oldroyd-B fluid in a half-space is readily es¬ 
tablished from the solution to the transient version of Stokes’ 
first problem for an Oldroyd-B fluid in a half-space, which is 
given by Eq. (7) in [13] (keeping in mind that a and ot, therein 
correspond to Xi and X 2 herein) using the convolution theo¬ 
rem for the Laplace transform. The resulting expressions are 
lengthy, hence, for the sake of brevity, we do not provide them 
here. (It is important to note, however, that the solution to the 
transient version of Stokes’ second problem for an Oldroyd-B 
fluid in a half-space given by Eq. (3.10) in [15] is erroneous 
for the reasons listed in [10]. Hence, carrying out the tedious 
integration and simplification may have scientific value.) 

On the other hand, noting that small-/ asymptotics corre¬ 
spond to large-p asymptotics in the Laplace transform domain 
(see, e.g.. Sec. 2.2 in [14]), we can obtain the early time be¬ 
havior of the transient version of Stokes’ second problem for 
an Oldroyd-B fluid in a half-space. To this end, specializing to 
the (5=0 case of Eq. (la) in [13], taking u(0,t) = H{t)cost 
as the first boundary condition in Eq. (lb) of [13] and applying 
the Laplace transform in /, we obtain the subsidiary ordinary 
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differential equation 


where % = a/Ai/A 2 and A = (Ai — X' 2 )l{X\Xi), whence 


( 1 ir. \ /A 1 ^ 

(TTI5j'’" = 37’ ’ 

»(Q^P) = 2 ^, I ’ «(«,/’) =0. (A.lb) 

/7‘^+ 1 

where i7 denotes the (temporal) Laplace transform of u with 
parameter p. Solving Eqs. (A.l) is straighforward: 





xexp{-yXv^} 



yjA 

2 VE 



(A.4) 


Inverting the Laplace transform term-by-term, we obtain 


u{y,p) = 


p ^+1 


exp < -jA 


' {1+Xip)p 

I+X 2 P 


Expanding for p ^ we find 


(A.2) 






-yjerfc 



(A.5) 


u{y,p) 



X exp 


-yxVp+ 


yx^ 

2 v^ 



as the early time (small-f) asymptotics. Clearly, the initial de¬ 
velopment of the boundary layer looks diffusive with a correc¬ 
tion proportional to A due to the non-Newtonian nature of the 
(p—>^ 00 ), (A.3) Oldroyd-B fluid. 
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